‘ — Z> < -3 


C T.t> <_~3 C=T~> <33 <3=3 <3=> C=3 <=3 C__> C=3 5. 

H M M hi M c=3 c=3 c,3 c_3 <r=> < 

A l_ ARSPrint 110972 

C:J ‘ J ‘ ^ [j\J \y J \s J j 

,'J= W M ’•'' | - 1 

u\J V/\J Cu J'J t/\i OflJ 


’ H '-- 


, .■=■> <. :=s <--.;> <_^> <„=> 

\J C/\J 0 v • 

M j» * f.^ij f.,4 :„*] « .& rj <u , 

C/V. v v"\j v/\/ v/v V/\i 




J A 


0 f 


'<(■ < /Of , .., ( ■ 

l/v V\j ./V VU oj VA 


AJ 


AJ -A; v 

j n- >■ r ^ ;■. 

/'J W Vb -.J v 

3 .' ~? <~‘-> r > c_ 




a!; CC; 


/j 


» c=r> c 33 > <33 


V, • / A A A 

‘ . V .. '3 f.V V'/ vrA 

• A A A A 

V -V . J O V <C> 

3. /' * 

' -••-'' An " Ax, 


> < . -> < 3 * 3 = 


> C= 3 > <33 • 


33> < — 3 c_3 ■ 


<33 <'33 <=> ■ 



> <3 3 C 3 ■ 


C 1.3 < 33 > CT 3 < 


t = 3 <. .3 1 


<3.3 C-. 3 <3. 3 «. ; 


ssssssssss 


.3 <-.-=> C.3 C= 



3 C 3 C.~3J c 3 ' 

/r- 


ssssssssss 
ssssssssss s s s s s s - 
ssssssssss sssssss- 
s I s s s s s s s s s s s s s s s s s 
ssssssssss 

S S S S S S S S S S S S S S S S S 2 

sss ssssss 



Application of Multispectral 

Remote Sensing to Soil 
Survey Research in Indiana 


■:Sr#;h 

>'.'/ 'J' y *J\i w A 

C A r,j ,7 t. 
3 - ("'ll '~'i 

W 3 -J bu 

' ' A'J Uj v\i 

A- A " s ? 

■■ A-. : /N . 

•. / ■■ • \. Kj\i 

,; ^TAiAA;^ 

A -'A v -.j j\ \j\'- 


s 

(T 0* (? 

An a 

c 2 ) qA oa <=_-a c-. 


o > ^ 'A/ O' ^ ^ ^ cyb Cr 'J ^ A •' l y 0 _ • |_ , . 

SS SSSSSSSSSSSSSS--x;A-jkxZachary,, J.JE.Cipra, 

R, I, Diderickson, S. J. Knstof 


: < .35 c. .. y ; ' 


SSS SSSSSSSSS SSSSSS Baumgardner 

S SS SSSSSS S MS 

SSSSSSSSSSS 


^ <33 CJ C.33 


t;<s « SSSSSSS— -SS 1 -*' 


! ' "I I r"> i -,'1 iA ' - ’ 

SS S 3 u-J vv/ 

n; 


a 

S A $ 

<3) 


A A A ^ c.-.- 3 


Sv U'u 1/<J 

13 v 


C. 3 
vA/ 


h2) A ^ 

sssssssssssss* 



S S S S - s s s s s s s 

C A A CC <c cs c St A 


, S3 ru s 

dj rJJ \AJ 3 

^ ; a a A A A A A - - A « - M r~ -j r,c , 

A cv )> ^ <-_A ’sA A-3) tr/ <s2t <s„A \A i v'V '/A 

f \ ^ A A ^ \A t A , C 

A;;oAji;s i 

(t ce 



C5 < . 5 < -5 1 


A A ^ Q: 


■C’o 

h) h) h) 's.A =s 2 / <iA A i/A l/\| 

S ^ ^ ^ ^ ^ ^ ^ ^ -Hr f, 'V! 

Ab C.A <SrV rb Ate ^ ^ S:A -^V ^ A UV UV iv/Vj 

\ A C C ^ vC h CC ^ C - 1 *uA A/ f n /) 

„A 'cb) Aa < 2 ) *<fc-S) v^) c 3 ) cA A'-./ L/vj vb 

^ ^ c <s c, <c S c. 3 h f'h 

■^33 Q 3 ^ <33 <S 3 > 03 / 3 A V'J/ lii/ -3 3 l/\; xtv 

(A* (A A’’ /(•' rr 3 G* fc* m* Sr* Sr* 

AA <iA o.A *t.A 


S S S A A A 
A A A A S A 

C;A 03 / < 3 / OA OA oA' 

CC CC CC 


C.C3 C„=> C 


'i=A 03) 

A 

oA 


CT ..3 C =3 < 3 =P 


S. SASSS^ 
SS A 


O «=_=> C3’ 


Purdue University 

<=33 C— > 1.3 > X > f_Z> < .><„><• .0 <.. -> ; _ _> * _. .3 : 

West Lafayette . Indiana 


c:> Ci c3.) e 


;; J; £ 


A - * r 1 s a n r. ( 

' v - ’/'-I 1 A AvA 

'b ’y,; s/\j uiJ 

■■ r 2 s :’ 1 f rj y 



Application of Mul tispectral Remote Sensing 
to Soil Survey Research in Indiana 1 

by 

A. L. Zachary, J. E. Cipra, R. I. Diderickson, 

S. J. Kristof, and M. F. Baumgardner 

ABSTRACT 

Recent advances have been made in the technology of measuring 
radiance from the earth's surface using multiple -walength airborne 
scanning spectrometers. Concurrently, advances were being made 
in the application of computer- implemented pattern recognition 
techniques to these multispectral data. 2 Together these two tools 
have resulted in a capability for mapping various earth surface 
features with extreme rapidity and varying degrees of accuracy. 

This study compared computer- implemented mappings based on spectral 
properties of bare soil surfaces with manning units of interest to 
soil surveyors. Some soil types could be differentiated by their 
spectral properties. In other cases, soils with similar surface 
colors and textures could not be distinguished spectrally. The 
spectral maps seemed useful for delineating boundaries between 
soils in many cases. 

Bowers and Hanks (2) measured laboratory reflectance in the 
400 to 2500 nm wavelength region of four Kansas soils. They con- 
cluded that surface moisture and organic matter strongly influence 
the reflectance and absorbance of solar radiant energy by soils. 
Cipra et al. (3) measured the reflectance of samples representing 
seven Indiana soil series under field conditions. They attributed 
percentage of visible incident energy reflected to soil color, 
texture, organic matter content, moisture content, and surface 
condition. The magnitude of the influence of each of these factors 
and their interactions were not discussed. Condit (4) examined 
spectral properties of 160 surface soil samples collected at various 
locations across the United States. Measuring laboratory condi- 
tions, he concluded that the general shapes of the spectral curves 
obtained for these soils could be classified into three types. 

These three types of curves could be represented by the chernozem- 
type soils, the pedalfer- type silts, and the laterite- type soils. 


'This work was sponsored under NASA Grant NGL 15-005-112 in 
cooperation with Purdue University and the Laboratory for Applications 
of Remote Sensing (LARS), West Lafayette, Indiana. 

laboratory for Agricultural Remote Sensing. 1970. Remote Multi- 
spectral Sensing in Agriculture, Vol . 4 (Annual Report). Research 
Bulletin No. 873, Agricultural Experiment Station, Purdue 
University, Lafayette, Indiana. 
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Baumgardner et_ al . (1) , using spectral measurements collected 
from aircraft, reported that soil organic matter appeared to be 
a dominant factor affecting reflectance when organic matter content 
was greater than 2 %. Their study was conducted on a 25-ha test 
site in Indiana which included eight soil series, using spectral 
measurements in 12 wavelength bands ranging from 400 to 2600 nm. 

Kristof (5) and Kristof and Zachary (6) conducted soils 
studies using multivariate pattern recognition techniques and 
computer processing of multispectral data collected by an airborne 
scanning spectrometer. Kristof and Zachary (6) concluded that 
"mapping" of soil types using these computerized procedures was 
partially successful. 

In the present study we have attempted to determine hoi«7 
favorably the spectral maps produced by computer processing 
comoared with conventional soil survey maps. Additionally, in 
cases where good agreement \vas not obtained, we have attempted 
to determine why the discrepancies occurred. 

Materials and Methods 

The three areas studied were designated as Soil Test 
Area 3 (STA 3), Soil Test Area 4 (STA 4), and Soil Test Area 5 
(STA 5). Soil Test Area 3 is located along U.S, Highway 37 
in Morgan County, in south central Indiana. The soils in STA 3 
were developed in late Wisconsin glacial material, including 
till, outwash, and aeolian sands. They are Alfisols (Gray-Brown 
Podzolic) and Mollisols (Humic Gley and Alluvial soils). Topo- 
graphy is nearly level to rolling. 

Soil Test Area 4 and 5 are located in Tippecanoe County, in 
west central Indiana. Soils in STA 4 are within the region of 
the Alfisols and include some wet Mollisols. These soils were 
developed in 45 to 90 cm of silt over glacial till. Soil Test 
Area 5 is also within the Alfisol region but the surface horizons 
are somewhat darker and contain slightly more organic matter 
than soils of STA 4. The area includes some wet Mollisols. The 
soils in the northern half of STA 5 were developed in moderately 
deep silts (1 to 1 1/2 m) ; whereas those in the southern half 
were developed in glacial till with less than 40 cm of silt at 
the surface. Table 1 gives the classifications of soil series 
occurring in STA 3, STA 4, and STA 5. 

The three study areas were field mapped at medium intensity 
using conventional soil survey procedures, giving considerable 
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attention to detail. Base photos we re color for STA 3 and black 
and white for the other two areas. Aerial Tnultisnectral scanner 
data were collected by the University of Michigan’s C-47 aircraft 
on April 28, 1967, at 1100 hours at an altitude of 1200 m above 
terrain (STA 3) and on May 26, 1969, at 1200 hours at an altitude 
of 1200 m above terrain (STA 4 and STA 5). Table 2 gives channel 
number disignations and wavelength bands for data collected. 

The multispectral scanner data in analog form were digitized 
and then analyzed using LARSYS programs. 4 Soil Test Area 3 was 
analyzed using the unsupervised classifier NSCLAS (8), which uses 
a clustering algorithm to classify spectral data into the number 
of classes specified by the researcher. This classifier is unsup- 
ervised in the sense that "training” areas are not input by the 
researcher, that is, the researcher does not define the "training" 
classes from the ground observations. He specifies only the 
wavelength bands to be used, the rectangular area or areas to 
be classified, and the number of classes. The resulting classi- 
fication then, is based entirely on the soectral data, and can be 
evaluated in terms of ground observations if desired. 

The channels selected for the map shown in Fig. 1 were 
2, 4, 6, 8, 12 and 12 (Table 2). After some experimentation it 
was decided to use 13 classes for the analysis of STA 3 and some 
surrounding area, which resulted in 8 classes being mapped within 
the boundaries of STA 3. This was done because programs allowed 
only rectangular boundaries to be entered into analysis procedures. 
Thirteen classes gave the best separation into two classes -- green 
vegetation and nonvegetated soils. 

Areas 4 and 5 were analyzed using a supervised classification 
approach (7). All eleven wavelength bands were used (Table 2). 

This classifier, $CLASSIFY, uses a maximum likelihood algorithm 
in the decision-making process. In this specific case, the 
analyst defines classes on the basis of field-observed soil 
types, and the computer uses spectral data from these "training" 
areas to characterize each soil type. 

Results 

Figure 1 is a computer printout of STA 3. The northeast 
area of Princeton fine sandy loam is represented predominantly 
by (.) with some (-) and (=) intermixed. In the middle of the 
printout Princeton fine sandy loam is a uniform area represented 
predominantly (.). The areas of Princeton fine sandy loam on 
the southwest side of the farm are represented predominantly by 


4 LARSYS is a software package developed by LARS for handling and 
analysis of multispectral data in digital form. A more complete 
description is given in the reference cited in footnote 3. 
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(-) and (=) . Only one soil type of Princeton was mapped, however, 
reexamination in the field revealed the northeast area was 
predominantly loamy fine sand with some inclusions of fine sandy 
loam, while the southeast area was predominantly fine sandy loam 
with some inclusions of loamy fine sand. The two areas of Princeton 
differed in slope and organic matter content as well as texture. 

Most of the area mapped Ockley loam on the soils map is 
represented by (=) on the printout. However, the east part 
of the Ockley is represented by (+). The large area mapped Fox 
loam is represented by (*) and (I) on the printout. Field 
examination revealed that the area with symbol (*) contained 
more sand in the surface than area represented by (I) but both 
areas fell within the range of characteristics of Fox loam. 

The east side of the large area mapped Ross is represented 
largely by (*) and (+) with some (I) present. The area where 
(*) appears was later found to be an inclusion of Fox loam, 
which contained more sand in the surface horizon than the area 
mapped Ross. 

The west side of the area mapped Ross is represented largely 
by (0) . Field examination revealed no reason why the computer 
printout showed two different symbols (0 and I) in the Ross 
area. The computer printout showed (0) and (H) for the area 
mapped Rensselear and the west part of the area mapped Ross. 

The soils of the glacial till area in the northeast part of the 
map (Miami and Crosby soils) were not well differentiated from 
the outwash and aeolian soils. 

Upon further examination of the spectral data it was found 
that the measured values of spectral response were less relaible 
in the right hand one-third to the computer map. This was 
because of sun angle and/or look angle effects causing an 
apparent darkening of the data. This effect is most apparent 
to the right of a line from the legend symbols "Ro" to the 
legend symbol ”F M . It can be noted there is little agreement 
between soil boundaries and computer mapped boundaries beyond 
this point. 

Figure 2 shows the soil man and the computer printout map 
for STA 4. 

The Russell soils of the computer printout, represented by 
(-), compared well with the Russell soils delineated on the 
soils map. The printout shows the delineations of the Russell 
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soils vary slightly from the soils map. In recheclcing these 
areas it was observed that the printout for the Russell series 
was more accurate than the soils map. 

Sandy areas in STA 4 were readily separated from silt loam 
and silty clay loam areas by the pattern recognition techniques. 

The area in the northeast corner represented by (/), the symbol 
for the fletea soils, correspond very well with the Metea delinea- 
tion on the soils map. Some other areas not mapped as Metea 
sandy loam on the soil map were indicated as Metea on the printout. 
Reexamination of these areas showed there was more sand in the 
surface 10-20 cm than in soils of surrounding areas. However, 
these areas are not classified as Metea soils, but are inclusions 
of other soils of STA 4. 

The area mapped Kokomo is rather uniform, except in some 
areas where some 1 ight- colored overburden has been mixed with 
the plow layer of the dark-colored Kokomo soil. In general, 
there was good agreement between the map and the printout for 
the Kokomo area. 

In the south part of the farm the area mapped Toronto was 
also well separated on the printout. However, on the west 
side of the ditch the printout indicated Toronto soil where 
Del Rey and Kokomo were mapped. The Del Rey soil west of 
the ditch has a darker surface than that which is described as 
modal for the series. The Kokomo surface horizon is lighter 
colored than modal for the series because of some mixture of 
light- colored soils deposited from higher topographic positions. 
This area probably showed up as Toronto on the printout because 
the surface properties of the Kokomo and Del Rey soils in this 
area are similar to those of Toronto. 

At the time this area was flown in May, part of the area 
west of the ditch was covered by oats approximately 20 cm tall. 

The computer was "trained” on separate samples in this area 
and was able to differentiate among series to some extent. It 
is not known at this time how much vegetative ground cover can 
be present without obscuring soil patterns. 

The soils of the STA 5 (Fig. 3) are mainly Ragsdale silty 
clay loam (Typic Argiaquoll) and Reeseville (Acric Ochraqualf) 
in the northern part. In the southern part soils are mostly 
Brookston silty clay loam and silt loam (Typic Argiaquoll), 

Crosby silt loam (Aerie Ochraqualf) , Celina silt loam (Aquic 
Hapludalf) , and Reeseville silt loam. 
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There is excellent agreement among some of the areas on the 
soils map and the printout (for example, the Reeseville in the 
northern part). The Ragsdale soil was also well-delineated on 
the printout except for some inclusions on Brookston silt loam 
and silty clay loam. These inclusions were also observed in the 
field, but were too small and intermixed to delineate on the soil 
map. In the southwest part of the test area there was good 
identification of Brookston silty clay loam on the printout. The 
Brookston silty clay loam area showed a small percentage of 
Brookston silt loam and Ragsdale silty clay loam on the printout. 
Field check verified some inclusions of these two soils. 

Reeseville soil in the southern part was accurately 
identified in the western part of the area, however, in the 
eastern part much of the Reeseville area was incorrectly 
identified by the computer as Crosby. Slightly sandier surface 
texture and darker color in the eastern part may have caused 
this problem. In this area the surface color of Crosby ranges 
from dark gray to grayish brown. Celina was delineated on the 
printout very well. Toronto was not well-delineated on the 
printout. In field mapping the areas of Toronto were small and 
hard to separate from Brookston and Ragsdale. Since Toronto 
is a transitional soil, it was difficult to distinguish from 
the darker-colored Mollisols by spectral properties and pattern 
recognition techniques. 

For further evaluation of multispectral remote sensing 
technology in soil survey, training samples were taken from 
STA 5 and and an attempt was made to extend the same mapping 
units beyond STA 5 using the computer. This was done for an 
area south of STA 5. The training samples were adequate for 
about 3 km. Beyond this distance discrepancies were noted between 
ground observations and computer identification of soils. 

Conclusions 

This study revealed a definite relationship between 
multispectral imagery and soil types. Supervised classifications 
gave results which agreed more closely with the soil survey map 
than did unsupervised classifications. Sun-angle or look-angle 
effects, or both, were believed to limit the sensitivity of the 
method. This effect is ussally less pronounced in data collected 
near solar noon. In spite of these present limitations, it is 
believed that multispectral remote sensing and computerized 
pattern recognition techniques have potential in the area of 
soil mapping. Large areas of bare soil can be "mapped” rapidly 
by computer techniques and these maps may provide the soil 
scientist a useful supplement to aerial photography when making 
soil surveys. 
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Table 1. Soils Occurring in the Study Areas. 


Soil Type 


Classification 


Crosby loam 
Miami loam 
Ockley loam 

Princeton fine sandy loam 
Martinsville loam 
Fox loam 
Ross silt loam 
Rensselear fine sandy loam 
Kokomo silty clay loam 
Brooks ton silty clay loam 
Toronto silt loam 
Mete a silt loam 
Del Rey silt loam 
Fincastle silt loam 
Xenia silt loam 
Russell silt loam 
Ragsdale silty clay loam 
Brookston silt loam 
Crosby silt loam 
Celina silt loam 
Reeseville silt loam 


Aerie Ochraqualf 
Typic Hapludalf 
Typic Hapludalf 
Typic Hapludalf 
Typic Hapludalf 
Typic Hapuldalf 
Cumulic Hapludoll 
Cumulic Hapludoll 
Typic Argiaquoll 
Typic Argiaquoll 
Udollic Ochraqualf 
Arenic Hapludalf 
Aerie Ochraqualf 
Aerie Ochraqualf 
Aquic Hapludalf 
Typic Hapludalf 
Typic Argiaquoll 
Typic Argiaquoll 
Aerie Ochraquoll 
Aquic Hapludalf 
Aerie Ochraqualf 



Table 2. 


Wavelength bands used in this study. 


STA 3 

Channel no. 
2 
4 
6 

8 


10 

12 


wavelength band (ra) 
440-460 
480-500 
520-550 
580-620 
660-720 
800-1000 


STA 4 and STA 5 
Channel no. 

1 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


wavelength band (nra) 
400-440 
520-550 
550-580 
580-620 
620-660 
660-720 
720-800 
800-1000 
1000-1400 
1500-1800 
2000-2600 




SOIL MAPPING UNITS 

C Crosby I Ma Martinsville I 

Mi Miami I F Fox I 

0 Ockley I Ro Ross sil 

P Princeton fsl Re Rensselear fsl 


Figure 1. Nonsupervised computer classification for STA 2 with 
soil survey map overlaid. Blank areas are primarily 
vegetation and other non-soil materials. 
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TRAINING CLASS SYMBOLS 

Z Fincastle sil 
Xenia sil 
- Russell sil 
I Toronto sil 
H Brookston sil 
H Brookston sicl 
M Kokomo sicl 
/ Metea sl 
= Del Ray sil 


SOIL MAPPING UNITS 

f:-: Fincastle sil 

EHZZZ3 Xenia sil 
fgglU Russell sil 
C2HSZ-) Toronto si I 
1 I Brookston si cl 

r l Kokomo si cl 

i Metea sl 
EHS3 Del Roy sil 


Figure 2 


Computer classification and soil survey map of STA 4. 
Blank areas on the classification indicate "threshold 
points; no classification decision made. 






